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Label 


Size 


$g 


P 




m uni. (MeV) 


# configs. 


m P n Q > (MeV) 


32ID 


32 3 x 64 x 32 


Iwasaki+DSDR 


1.75 


2 


170, 250 


181, 148 


140 


241 


24 3 x 64 x 16 


Iwasaki 


2.13 


2 


330, 420 


202, 178 


240 


321 


32 3 x 64 x 16 


Iwasaki 


2.25 


3 


290, 350, 400 


300,312, 252 


220 



Table 1: A summary of the properties of the three ensemble sets used in this analysis. Here 'Sg' denotes the 
gauge action, W ens ' the number of ensembles, 'm^. nL ' the unitary pion mass on each of those ensembles, '# 
configs' the number of gauge configurations used in this analysis, and m P K Q > the lightest available partially- 
quenched pion mass on the ensemble set. 

1. Introduction 

In these proceedings we present continuum limit predictions for the neutral kaon mixing pa- 
rameter B K , the strange and average up/down quark masses and the pion decay constant f % , ob- 
tained via simultaneous chiral/continuum fits to the three domain wall fermion ensemble sets given 
in table [I]. Our most recently generated ensemble set, referred to by the label '32ID', uses a mod- 
ified Iwasaki gauge action that includes the Dislocation Supressing Determinant Ratio (DSDR) 
term, which allows us to simulate with near-physical pion masses on a coarser j8 = 1.75 lattice, 
while retaining good chiral symmetry and topology tunneling. 

Combined analyses of the 241 and 321 lattices, which we refer to collectively as the '2010 anal- 
ysis', have recently been published §]. In these papers we developed a strategy for the combined 
analysis of multiple ensemble sets that maximises the use of the available data in constraining the 
fits. In these proceedings we develop this strategy further to include 32ID ensemble set. Note that 
the number of configurations available on the lightest 32ID ensemble set has almost doubled since 
the conference, and the heavier ensemble has also increased in size by 35%, hence the conclusions 
presented here differ slightly from those given at the conference. 

The layout of these proceedings is as follows: We provide a short discussion on the DSDR 
term, followed by a summary of our combined fitting strategy. We then present results for the pion 
decay constant before briefly touching on the non-perturbative renormalisation techniques used in 
the calculations of B K and the physical quark masses. Results for these quantities follow. We 
conclude with a brief summary and outlook for this analysis. 

2. The DSDR Term 

The 2010 analysis was performed to data over the range m K = 290 — 420 MeV. Here the extrap- 
olations down to the physical pion mass of ~ 135 MeV provided the dominant contribution to the 
systematic errors on the continuum predictions. For example, the chiral extrapolation systematic 
on f K was ^ 4%. This provided a strong motivation for reaching down to lighter quark masses. 

In order to avoid finite- volume effects, simulations with lighter quark masses require lattices 
with a larger physical volume. However, at our typical couplings the required increase in the 
number of lattice sites is beyond the reach of our current computing resources; we are therefore 
forced to simulate with coarser lattices. This has the unfortunate side-effect of increasing the size of 
the chiral symmetry breaking effects in the domain wall fermion formulation, due to the following 
mechanism. 

The size of the chiral symmetry breaking is parameterised by an additive mass-shift known 
as the 'residual mass' m res . This quantity is governed by the eigenvalue density p of the four- 
dimensional Hamiltonian H? = 2tanh _1 ( 2 +l!v ) describing quark propagation through the fifth 

dimension. Here D w is the Wilson Dirac operator and H w = ^D w is the hermitian Wilson Dirac 
operator. The equation above implies a relationship between the eigenmodes of H T and H w . The 
modes of the latter are divided into two regions by a mobility edge X c \ those above X c are extended 
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h — h e=0.040, ^=0. 
q-0 e=0.005, e =0. 
H-g e=0.010, £,,=0.100 




Figure 1: An illustration of the molecular dynamics force imparted by the DSDR term at various values of 
£f and £t>, reproduced from ref. ^ 

over the lattice whereas those below are localized - in fact this property is essential in ensuring 
that the four-dimensional effective theory of the fields on the boundary is local. The structure of 
eigenmodes of H w implies that m res has the following dependence on L s [0]: 

m res ~ R 4 eP (X c ) t^L + Rjp (0) 1 (2. 1) 

where the first term contains contributions from the extended modes and the second term from the 
localized near-zero modes. 

In modern simulations, L s is typically large enough that m res is dominated by the near-zero 
mode contribution. These modes are associated with small tears or 'dislocations' in the gauge 
field, which occur more often as we approach the disordered strong-coupling region. In order to 
retain good chiral symmetry at stronger coupling we must therefore seek to suppress these modes. 
However we must retain enough of the very-near-zero modes that allow topological tunneling to 
occur. This can be achieved by introducing a weighting term, the DSDR term, into the gauge 
action, given by ^ g |] 

det[D w (-M + /8^)tp w (-M + /8^)] _ rr ^ + 8 / 2 
1 ' f ' b) ~ det[D w (-M + ie f fYD w (-M + i£ f f)} "ll^ + e, 2 ' ^ 

where A; are eigenvalues of H w and £f and £^ are tunable parameters. This introduces a force in 
the molecular dynamics evolution of the form 



&i(e f ,e b ) = — -log-L— £ , (2.3) 




which can be tuned to peak in the near-zero region without further suppressing the very-near-zero 
modes. In figure [l] we reproduce a plot from ref. [[|] which shows the force as a function of X for 
several combinations of £f and on a test simulation, demonstrating this suppression. 

3. Simultaneous Fitting Procedure 

We obtain our fit forms via a dual expansion to next-to-leading order (NLO) in a 2 and the light 
and heavy quark masses m/ and nth, adopting a power counting that discards terms like &(a 2 m) and 
&{m 2 ) and higher at this order. Expanding around a non-zero mass point (m/o^/^o) and absorbing 
constant terms containing m/o into the leading coefficient, we obtain the analytic fit function 

/// =cf'(l + C f a 2 ) + Cf ^±Hhl + c^fht + C 3 A (m h - m h0 ) (3.1) 
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2010 analysis 


This analysis (all data) 


This analysis (m n < 350 MeV) 


fn (MeV) 


124(2)(5)(2) 


125(2)(2)(2) 


127(3)(0.5)(3) 


m^S(2 GeV) (MeV) 


3.59(13)(12)(6)(8) 


3.48(6)(7)(3)(8) 


3.39(9)(4)(2)(7) 


mfS(2 GeV) (MeV) 


96.2(1.5)(0.2)(0.1)(2.1) 


94.9(1.2)(1.4)(0.2)(2.1) 


94.2(1.9)(1.0)(0.4)(2.1) 


B K 


0.749(7)(21)(3)(15) 


0.748(6)(15)(4)(15) 


0.751(11)(8)(4)(14) 



Table 2: Results for f K , the average up/down quark mass, the strange quark mass and Bk- The first column 
contains the 2010 analysis result [J, the second the result obtained in this analysis by fitting to the full 
range of available data, and the third by fitting only to data with m K < 350 MeV The errors are statistical, 
chiral, finite- volume and NPR errors (where appropriate) respectively. 



for the pion decay constant and similar forms for other quantities. Here m = m + m res . We also 
consider expanding about the SU (2) chiral limit (with and without finite-volume corrections), in 
which case we obtain the usual NLO ChPT fit forms with an extra a 2 and (nth — mno) term. We 
refer to these as the ChPT and ChPTFV fit forms. Using these three fit ansatze, we simultaneously 
fit the 321, 241 and 32ID ensemble sets, with the coefficients shared between all three data sets. 
For ensemble sets i other than the primary set, chosen as the 321 set, we must include appropriate 
factors of the ratio of lattice spacings R a = a 321 /a 1 and quark masses R t / h = m^jl/m)^ in the fit 
forms such that the coefficients C can be assumed equal. In this analysis we adopt the so-called 
'generic scaling' approach ^ in which these ratios are determined as free parameters in the fit. 

Before taking the continuum limit, we must fix the primary lattice spacing a 321 and the quark 
masses m 32l d and m 321 in 321 normalisation. This is achieved by varying these parameters until the 
continuum values of the pion, kaon and Omega baryon masses agree with their physical values. 
These quantities define the 'scaling trajectory' along which the continuum limit is defined. Note 
that this makes the a 2 coefficients of these quantities zero by definition. As the above conditions 
are applied in the continuum limit, the procedure is necessarily iterative. 

The inclusion of the Iwasaki+DSDR data complicates the situation slightly over the 2010 
analysis, in that the coefficients of the a 2 terms must now be allowed to differ between the two 
different gauge actions. Although this introduces one extra degree of freedom per quantity into the 
fits, any algorithmic instability that this may cause is offset by the increased number of data points. 

Recall that eqn. [O] contains a term in the heavy sea-quark mass m^. We vary this parameter 
using reweighting [Q, 0], whereby the weight of a given gauge configuration is re-evaluated in the 
path integral at several strange quark masses shifted by up to 20% from the simulated value. This 
allows us to explore the strange sea quark dependence with minimal cost and to ultimately quote 
predictions at the physical strange quark mass at the penalty of an increase in statistical error. 

4. Preliminary Results For f K 

f K is obtained from the (0|Ao|ft) matrix elements in the usual way, the only difference being 
that for domain wall fermions one must correctly renormalise the 4d axial current to match the 
continuum current. Following ref. the renormalisation factor is obtained from the ratio of the 
5d DWF conserved vector current to the 4d vector current: this was shown to be more precise than 
the ratio of axial currents due to the unknown renormalisation coefficient between the 5d DWF 
PC AC current and the continuum current. 

Figure || shows the chiral extrapolation of f n down to the physical up/down quark mass using 
the analytic and ChPTFV ansatze. Following the 2010 analysis we estimate the error on the chiral 
extrapolation through the difference of the ChPTFV and analytic predictions, and the finite-volume 
error from the difference of the ChPTFV and ChPT predictions. We obtain the value given in the 
second column of table ||. This result is some 4% (1.5a) below the physical value of 130.7 MeV. 
A similar discrepancy was noted in the 2010 analysis (first column of the table), and was attributed 
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mi (GeV) mi (GeV) 



Figure 2: Comparisons of the unitary f n data corrected to the continuum limit and fit using the ChPTFV and 
analytic ansatze, where the fits are performed to the full data set (left) and only to the data with m K < 350 
MeV (right). The data included in the fit are marked with circles, and those excluded with diamonds. The 
continuum limits are marked with squares. Points in pastel colours are corrected to the continuum using the 
ChPTFV ansatz, and in bold colours by the analytic ansatz. 

to the systematic error on the chiral extrapolation. There, the strategy of estimating the chiral 
error from the difference of the ChPTFV and analytic results produced a systematic error sufficient 
to explain the discrepancy, but for this analysis it appears to be an underestimate. However, the 
introduction of the light DSDR data allows us to perform stable fits even after removing some 
of the heavier data from the 321 and 241 ensemble sets. We can therefore restrict our fits to a 
region of lighter mass in which we would expect the fit ansatze to perform better. Cutting the 
heaviest two ensembles (cf. table []]) such that the heaviest pion has a mass of 350 MeV, we obtain 
the result given in the third column of table [| Here, as a result of the cut, the central value has 
increased such that the result is now consistent with the physical value, even without incorporating 
the chiral systematic. Notice also that our ad hoc chiral systematic has decreased substantially, and 
the # 2 /d.o.f. also decreased, which suggests that the fits do indeed perform better when restricted 
to this lighter mass range. Although this analysis is still preliminary, these results are promising. 

5. Renormalisation for Bk and the quark masses 

Before giving results for the quark masses and B K , we present a short summary of the non- 
perturbative renormalisation techniques that are used to convert these quantities into the canonical 
MS scheme with unprecendented precision. 

Direct conversion into the MS scheme on the lattice is not possible as this scheme is regu- 
larised in non-integer dimensions. Instead we first convert to a convenient intermediate scheme, 
which is then run to high energy and matched to MS using perturbation theory. We use variants 
of the Rome-Southampton Regularisation-Invariant Momentum (RI-MOM) scheme, in which the 
renormalisation coefficients are defined from the ratio of the amputated Green's function of an op- 
erator to its tree-level value in the Landau gauge at a particular momentum scale. In the RI-MOM 
scheme, the Green's functions are formed using massless propagators of equal momentum p, and 
the scale is defined as ji 2 = p 2 . However this is a so-called 'exceptional' momentum configuration 
in which the hard external momenta can be routed in such a way that parts of the diagram contain 
only soft momenta - this enhances the effect of the spontaneous chiral symmetry breaking at high 
momenta. In order to avoid this problem we follow the 2010 analysis in adopting 'symmetric' 
momentum conditions (SMOM), for which the external momenta are not equal but rather obey 
/I 2 = Pi= p 2 = (pi — Pi) 1 - Note that this applies both for the mass renormalisation, which we 
obtain from the scalar bilinear vertex, and the four-point VV +AA weak operator contained in the 
definition of B K ; for the latter the Wick contractions comprise traces of two bilinear vertices - the 
momenta p\ and p2 are assigned to the two propagators forming each of those bilinears. 
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As in the 2010 analysis, we use volume source propagators in these calculations, which offer 
a significant improvement over the traditional point sources, giving statistical errors of the order 
0.1% even with only <^(20) gauge configurations. 

One of the dominant systematic errors on the NPR in our earlier analyses arose through the 
use of momenta whose unit vectors were not equivalent under hypercubic rotations, and hence do 
not have the same discretisation errors; this induces a scatter in the RI-SMOM renormalisation 
coefficients as a function of momentum. In ref. ^ we described how this can be corrected using 
twisted boundary conditions in the valence sector, which allow us to smoothly vary the magnitude 
of the momentum while keeping the direction fixed. We also discussed how another of the dominant 
systematic errors, that associated with the truncation of the perturbative series, can be reduced by 
a factor of two by performing the matching to the MS scheme at 3 GeV rather than the canonical 2 
GeV. For this analysis these improvements have been applied in the case of Bk but have yet to be 
applied to the quark mass renormalisation. 

6. Results for the Quark Masses 

The continuum physical quark masses are obtained in the normalisation of the 321 ensemble 
set, and hence must be renormalised into the MS scheme to remove any cutoff dependence. The 
renormalisation coefficients are determined by performing the continuum extrapolation of Z m /Ri/ h 
over the two Iwasaki lattices, where R t / h are the quark mass ratios defined in section || and R 3 ^ = 1 
by definition. Notice that the coefficient on the 32ID lattice is not needed for this procedure. As 
mentioned above, the MS (3 GeV) lattice coefficients with twisted boundary conditions have yet 
to be determined, hence for this analysis we reuse those given in ref. [||], although the continuum 
extrapolation is performed anew with the lattice spacings and R^ h obtained here. Following the 
2010 analysis procedure, we choose our best NPR scheme for the central value take the truncation 
error on the renormalisation from the size of the two-loop contribution to the MS matching. 

From a fit to the full data set, we obtain the values given in the second column of table ^. 
For comparison, we give the 2010 analysis result in the first column. We find results that are very 
consistent, and observe a factor of two reduction in the statistical and chiral systematic errors on 
the up/down quark mass over the 2010 analysis as a result of including the lighter data. We also 
see a reduction in the statistical error on the strange quark mass, but also a large increase in the 
chiral error. This is likely a result of allowing the mass ratios Z/ and Z^ to differ between the fit 
ansatze, where before they were fixed to values obtained by matching the lattices at an unphysical 
mass scale (the fixed trajectory method). We intend to investigate this further. 

We also investigate the effect of cutting out the heaviest two ensembles on the quark masses. 
We obtain the result given in the third column of the table. Here as with f n , we see significant 
improvements in the estimated chiral error, at the expense of an increase in statistical error. 

7. Results for B K 

We obtain results for Bk through an independent simultaneous fit over our three ensemble 
sets, using the lattice spacings, quark masses, etc., obtained from the fits above. In particular, we 
constrain the ChPT fits by including the lowest-order ChPT parameters / and B from the main 
analysis. The fits are performed to MS-renormalised data, where we use the coefficients obtained 
by matching to our lattice scheme at 3 GeV as discussed above, after which we convert to the RGI- 
scheme for the convenience of the reader. The NPR error is obtained by taking the difference of 
our two best schemes, as discussed in ref. From the fits to the full data set, we obtain the result 
given in the second column of table ||. This result is very consistent with the 2010 result, and shows 
a 30% improvement in the chiral error. Fitting only to the data with m K < 350 MeV, we obtain the 
result given in the third column of the table. As before we see a substantial decrease in the chiral 
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mi (GeV) mi (GeV) 

Figure 3: Comparisons of the unitary Bk data in the SMOM(^,^) scheme corrected to the continuum limit 
and fit using the ChPTFV and analytic ansatze, where the fits are performed to the full data set (left) and 
only to the data with m n < 350 MeV (right). The symbols and their colours are described in the caption of 
figure |. 

error, but here the increase in the statistical error is larger than before. This is likely due to the lack 
of statistical resolution on the lightest 32ID data points. 

8. Conclusions and Outlook 

Using the Iwasaki+DSDR gauge action, we have been able to simulate with near physical 
pions while retaining good topological tunneling properties and small finite-volume corrections. 
Including these data in simultaneous fits with our Iwasaki lattices, we have been able to substan- 
tially improve our continuum predictions over the 2010 analysis, especially after cutting out the 
heaviest two ensembles such that the heaviest pion now has a mass of only 350 MeV - performing 
this cut we observed factor of two reductions in our estimated chiral systematic. In addition, our 
preliminary prediction for f % , 127(4) MeV, is now consistent with the physical value. We observe 
that NLO chiral perturbation theory extrapolations over the mass range of 140 — 350 MeV should 
be expected to fail at the 5% level, hence we expect this value to rise further towards the physical 
value as we further restrict the fit range in future analyses. 

It is our intention to publish an analysis of these data shortly, after which, alongside continuing 
to generate more data on the 32ID ensemble set, we intend to commence the generation of further 
domain wall fermion ensembles with near physical pions, taking advantage of the large increase in 
computing power provided by our upcoming IBM Blue Gene/Q resources. This will allow us to 
further refine our continuum predictions. 
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